Curvature corrections and topology change transition 
in brane-black hole systems: A perturbative approach 
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We consider curvature corrections to static, axisymmetric Dirac-Nambu-Goto membranes em- 
bedded into a spherically symmetric black hole spacetime with arbitrary number of dimensions. 
Since the next to leading order corrections in the effective brane action are quadratic in the brane 
thickness I, we adopt a linear perturbation approach in £ 2 . The perturbations are general in the 
sense that they are not restricted to the Rindler zone nor to the near-critical solutions of the un- 
perturbed system. As a result, an unexpected asymmetry in the perturbed system is found. In 
configurations, where the brane does not cross the black hole horizon, the perturbative approach 
does not lead to regular solutions if the number of the brane's spacetime dimensions D > 3. This 
condition, however, does not hold for the horizon crossing solutions. Consequently we argue that 
the presented perturbative approach breaks down for subcritical type solutions near the axis of the 
system for D > 3. Nevertheless, we can discuss topology-changing phase transitions in cases when 
D = 2 or 3, i.e. when the brane is a 1-dimensional string or a 2-dimensional sheet, respectively. 
For the general case, a different, non-perturbative approach should be sought. Based on the energy 
properties of those branes that are quasi-statically evolved from the equatorial configuration, we 
illustrate the results of the phase transition in the case of a D — 3 brane. It is found that small 
thickness perturbations do not modify the order of the transition, i.e. it remains first order just as 
in the case of vanishing thickness. 

PACS numbers: 04.70.Bw, 04.50.-h, 11.27.+d 



I. INTRODUCTION 

Higher dimensional black objects and branes are of im- 
portance and interest in several areas of present days 
physics. The classical black hole uniqueness theorems 
are known to fail in higher dimensions, and it turns out 
that a whole menagerie of black objects (black strings, 
rings, cigars, etc.) appear to exist. The study of new 
types of black objects became a very active research field 
recently, and among many other interesting aspects, the 
properties of possible transitions between the different 
types, or phases, is of special interest (see e.g. [l|, HJ and 
references therein). For example, during the transition 
between a caged black hole and a black string phase, Kol 
H demonstrated that the Euclidean topology of the sys- 
tem changes. This type of transition is called merger 
transition, and Kol found a strong similarity in its prop- 
erties with the Choptuik critical collapse phenomena [J] . 

Recently, Frolov suggested a simple toy model with 
many features in common with merger and topology 
changing transitions Q . The model consists of a bulk TV- 
dimensional black hole and a test 73-dimensional brane 
in it (D < N — 1), called brane-black hole (BBH) sys- 
tem. The black hole is spherically symmetric, static and 
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can be neutral or charged. The brane is infinitely thin, 
and it is described by the Dirac-Nambu-Goto p, 0, HI 
action. It is also static, spherically symmetric and it 
is assumed to reach asymptotic infinity in the form of 
a (D — l)-dimensional plane. Due to the gravitational 
attraction of the black hole, the brane is deformed and 
there are two types of equilibrium configurations. The 
brane either crosses the black hole horizon, or it lies to- 
tally outside of the black hole (see FIG. 1). In between 
the two types of configurations there exists a critical so- 
lution that separates the two phases. Frolov studied the 
transition between the so called subcritical phase (when 
the brane does not intersect the black hole horizon) and 
supercritical phase (when the brane crosses the horizon) , 
and found a close similarity both with the merger transi- 
tion in a caged black hole - black string system and with 
the Choptuik critical collapse phenomena. 

The AdS/CFT correspondence Q also provides moti- 
vation to study the above BBH system. In fact, according 
to the correspondence, at sufficiently high temperature, 
a small number of flavors (Nf) of fundamental matter in 
strongly coupled gauge theories with a large number of 
colors (N c ^S> Nf), may be described, in the holographic 
dual, by probe Dq-branes in the gravitational background 
of a black hole [lfl [U 0. In [lS GJ], using the tool of 
the gauge/gravity correspondence, Mateos et al. stud- 
ied the phase transition of quark-antiquark bound states 
(mesons). Their model was very similar to Frolov's toy 
model, and using the results of [f|, they demonstrated 
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that in the case of an infinitely thin brane, the system 
generally undergoes a first order phase transition char- 
acterized by a change in the meson spectrum. The cor- 
responding phase diagram in the vicinity of the critical 
solution exhibits a self similar structure, and this criti- 
cal behavior and the first order transition are essentially 
universal to all Dp/Dq systems. 

In [3] it was also pointed out, that higher order cor- 
rections to the brane effective action may cause, in prin- 
ciple, modifications to the above picture and it is likely 
that they spoil the system's scaling symmetry and self- 
similar behavior. Indeed, higher-derivative corrections to 
the D-brane action correspond to finite 't Hooft coupling 
corrections in the holographic dual, and provide a more 
realistic description of the system. These corrections may 
become important in the vicinity of the phase transition, 
since the curvature of the brane becomes large there. 

In the context of low-scale gravity theories, the possi- 
bility that a micro black hole may form in high energy 
collisions, like those at the LHC, re-creates a setup sim- 
ilar to the BBH system described above. In particular, 
the question whether a black hole may escape into the 
extra dimensional bulk has raised some attention, due to 
the potential phenomenological relevance [l5[ . Clarifying 
the role of the thickness of the brane in that context is 
also an important issue. 

The dynamics of branes keeps also attracting attention 
in the context of higher dimensional generalizations of 
the Bernstein conjecture QjJ and the study of the 
stability of brane-black hole systems [l8[ . 

For all the above reasons, it is important to go be- 
yond the approximation of zero thickness and consider 
higher order, curvature corrections coming from small 
thickness perturbations in the BBH system. Curvature 
corrections to the dynamics of domain walls without self- 
gravitation, in the case of non-zero thickness have been 
investigated earlier by Carter and Gregory [l9|,[2(|. They 
demonstrated that the next to leading order contribution 
is quadratic in the wall width (the brane thickness) and 
they obtained an exact, analytic expression for the cor- 
responding effective action in terms of the intrinsic Ricci 
scalar R and the extrinsic curvature scalar K . 

In a recent paper [2l[ , Frolov and Gorbonos studied the 
role of curvature corrections on topology changing tran- 
sitions based on the effective brane action presented in 
0, H3] ■ In this work the authors focused their attention 
to the near-critical solutions and similarly to 0, they 
restrict their investigations to the Rindlcr zone, i.e. very 
close to the black hole horizon, where the radius of the 
intersection of the brane with the bulk horizon is much 
smaller than the radius of the bulk horizon. As an in- 
teresting result they found that "the second order phase 
transition in such a system is modified and becomes first 
order". This, however, seems to be in contradiction with 
the results of Mateos et al. (see e.g. [II]), where the 
authors demonstrated that the phase transition in the 
unperturbed system is generally a first order one. Ad- 
ditionally, they find that when the spatial dimension of 



the brane is larger than 2, supercritical solutions behave 
quite differently from subcritical ones, and for supercriti- 
cal solutions there is no singularity resolution. According 
to their numerical analysis, they did not find evidence for 
the existence of such solutions. A possible explanation 
for this is that stiffness correction to the brane action 
break the symmetry between the super- and subcritical 
solutions 

In the present paper we re-consider thickness pertur- 
bations to the BBH system. We proceed within a more 
general framework than that of [21(. We consider the 
same curvature corrected effective brane action obtained 
by Carter and Gregory [2(| , but we do not restrict 
ourselves to work in the Rindler zone and in the near- 
critical solution region. In addition, we choose to follow 
a different path to obtain the dynamical equation for the 
perturbations. 

The plan of this paper is as follows. In Sec. II we 
make a quick overview of the model in the infinitely thin 
case, and reintroduce the BBH setup to make the paper 
self-contained. In Sec. Ill the curvature quantities are 
discussed, while in Sec. IV we obtain the Euler-Lagrange 
equation for the curvature corrected brane action. As for 
the latter, we follow the perturbative treatment by Carter 
and Gregory [2(| , in the sense that we treat the curvature 
corrections as small perturbations in the effective action. 
They are indeed very small as being quadratic in the per- 
turbation parameter. To obtain the dynamical equation 
for the perturbations, we use the quadratic perturbation 
parameter to expand the 4t/i-order Euler-Lagrange equa- 
tion and keep the linear terms only. As a result, a second 
order, linear equation is found to describe the perturba- 
tion function tp, with a very complicated source term. 
In Sec. V wc analyze the asymptotics of the perturba- 
tion equation, and find that there is no regular subcrit- 
ical solution on the axis of the system above a certain 
dimension. This implies that our perturbation method 
is not appropriate in this region and a non-perturbative 
solution may be in order to find the general solution of 
the problem. In sec. VI we write down the full Euler- 
Lagrange equation, but due to its very complicated and 
highly nonlinear form we do not discuss its solution in 
the present paper. In Sec. VII we present the analytic 
solution for far distances and the numerical solution in 
the near horizon region of the perturbation equation for 
various dimensions. In Sec. VIII we address the question 
of the phase transition in the case of a D = 3 dimensional 
brane. For this purpose we use the approach of Flachi et 
al. [22| based on the energy properties of a quasi-static 
brane evolution from the equatorial configuration. 



II. THE THIN BRANE MODEL 

Let us overview, in this section, the important proper- 
ties of the BBH system, introduced in [j| , that we intend 
to study in the presence of a small brane thickness in 
the following sections. We consider static brane config- 
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urations in the background of a static, spherically sym- 
metric bulk black hole. The metric of an TV-dimensional, 
spherically symmetric black hole spacctime is 

ds 2 = g ab dx a dx b = -fdt 2 + r x dr 2 + r 2 d£l 2 N _ 2 , (1) 

where / = f(r) and dil%_ 2 is the metric of an TV — 
2 dimensional unit sphere. One can define coordinates 
8i{i = 1, . . . , TV — 2) on this sphere with the relation 



dn 2 +1 = d9 2 +1 + sin 2 9 i+1 dQ 2 



(2) 



The explicit form of / is not important, it is only assumed 
that / is zero at the horizon tq, and it grows monotoni- 
cally to 1 at the spatial infinity r — ► oo, where it has the 
asymptotic form [23j, 



/=l-(^ 

r 



N-3 



(3) 



In the zero thickness case, the test brane configura- 
tions, in an external gravitational field, can be obtained 
by solving the equation of motion coming from the Dirac- 
Nambu-Goto action 0,0, H, 



d%/-det 7 _ , 



S 



where 7 is the induced metric on the brane 



lu,v = 9ab 



dx a dx b 
WW 



(4) 



(5) 



and ^(fi = 0, . . . , D — 1) are coordinates on the brane 
world sheet. The brane tension does not enter into the 
brane equations, thus for simplicity it can be put equal 
to 1. It is also assumed that the brane is static and 
spherically symmetric, and its surface is chosen to obey 
the equations 



'JV-2 



tt/2 



(6) 



With the above symmetry properties the brane world 
sheet can be defined by the function 9d-\ = 9(r) and we 
shall use coordinates £ M on the brane as 



C = {t 



> ' > Vli ■ 



6„} with n = D — 2 



(7) 



With this parametrization the induced metric on the 
brane is 



(8) 

where, and throughout this paper, a dot denotes the 
derivative with respect to r, and the action reduces 
to 



S = AtA n j C dr , 
Co = r n sin" 6\ll + frH 



(9) 
(10) 



where At is the interval of time and A n = 2n"/ 2 /T(n/2) 
is the surface area of a unit n-dimensional sphere. 

The brane configurations are the solutions of the Euler- 
Lagrange equation 



d_ i d Cp 
dr I d'O 



dC 

de 



= 



(11) 



B 3 3 + B 2 9 2 + Bi 9 + Bo =0 , (12) 



which for the Lagrangian (fT0|) reads 



Bo = - 



ncott 



Bi = —n cot ( 



Bi = 



B*=r 



r 

rf 



+ L 
/' 

+ (n + l)/ 



(13) 



For the supercritical case, a regular solution of ([12]) has 
the expansion near the horizon 



h + 9 (r-r ) + . 



with 8q 



ncot ( 



f> 



■ 2 



(14) 



hence it is uniquely determined by the initial condition 
Oo- In the subcritical case the brane does not cross the 
horizon, and its surface reaches the minimal distance 
from the black hole at r± > tq which, for symmetry rea- 
sons, occurs at 9 = 0. A regular solution of (fT2|) near this 
point has the behavior 

;r-n) 3/2 + ... , 



= T)y/r — n 



(15) 



where 



'2(n+ 1) 



B 3 (r) 



(16) 



48(n+3)r;r a / 



rj 2 r 



6?y 2 (n + l)r/ 2 



24r/ + / 



4(n(rr) 2 + 6) + 12) 



+ 3rV(2(». + 2)/ + r/) 



(17) 



and hence the solution is uniquely determined by the pa- 
rameter r\. In Q the coefficient a is not considered, 
however it is necessary to fix its value to get a regular 
solution of (Tj"2"|) , and it also appears later on in the per- 
turbation equation. 

The far distance solution of (fT2|) can be obtained from 
the condition that the brane behaves asymptotically as a 
(D — l)-dimensional plane. The solution can be searched 
in the form 



= 2 + ^( r ) ' 



and one can get for v{f) 
v{r) = 



E + £- if n >l, 



(18) 



(19) 
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action (1201) reduces to 




FIG. 1: The picture shows a sequence of subcritical and su- 
percritical thin brane solutions in the case when a D = 4 
dimensional brane embedded into alV = 5 dimensional bulk. 
The different configurations belong to different initial condi- 
tions of n and 9o, respectively. For simplicity, the bulk black 
hole's horizon radius is put to be 1, and R and Z are the 
standard cylindrical coordinates. 



where p can be referred as the distance of the brane from 
the equatorial plane at infinity, and both p and p' are well 
defined continuous functions of the initial parameters 9q 
or r\. 



III. CURVATURE CORRECTIONS 

Small thickness perturbations to the brane dynamics 
are derived from higher-order, curvature corrections to 
the effective action of the brane. In the present model 
we do not consider the self-gravitation of the brane, hence 
the curvature scalars are completely determined by the 
embedding black hole spacetime. In the approximation 
when all the relevant dynamical length scales L of the 
system are very large compared to the parameter i that 
characterizes the thickness of the brane, Carter and Gre- 
gory obtained the following exact, analytic expression for 
the effective action of the brane dynamics [2fJ , 



S = I d n (J-det % , 



s 4d 

M 1 



dR + C 2 K- 



, (20) 



where R is the Ricci scalar, K is the extrinsic curvature 
scalar and the coefficients C\ and C2 are expressed by 
the wall thickness parameter as 



Co = 



1 



24 3 
The parameter /1 is related to the thickness as 

1 



p,V2\ 



(21) 



(22) 



which originates from a field theoretical domain-wall 
model, where /i is the mass parameter and A is the cou- 
pling constant of the scalar field. 

Now let us consider the N dimensional BBH setup de- 
scribed in Sec. HH where the thickness corrected effective 



MAn £ dr , 



C = -^-Coll + eS] 



where we introduced the notations 



L 2 ' 



R 



C 2 L 2 
I 2 



K 



(23) 
(24) 

(25) 



In order to calculate the curvature scalars we follow the 
method described in [2(| HH using the Gauss-Codazzi 
formalism to split the quantities into their components 
orthogonal and parallel to the brane world sheet. 

Let n a be a unit geodesic vector field orthogonal to 
the brane and z the length parameter along the integral 
curves of n a . Then each constant z hypersurface has an 
intrinsic metric h a b and extrinsic curvature K a b defined 
by 



V.g ah V a zV 6 z 

Kb = 9ab - n a n b , 
Kab = h c a V c n b , 



(26) 

(27) 
(28) 



where the length parameter in our chosen spherical coor- 
dinate system is 



9ab^^dr^ 



f 



2 9' dr, (29) 



t is the curve parameter and the prime denotes the 
derivative with respect to r. We make the assumption 
that the perturbed system preserves the original sym- 
metry properties of the unperturbed BBH setup, hence 
expressing 2 as a function of r we find that the only iron 
vanishing components of n a are 



frO 



1 



n 



, - 2 



(30) 



l + fr 2 9 ryi + /r 2 6 

The extrinsic curvature scalar can be obtained as 

K ee K a a = h c a V c n a , (31) 

and the Ricci scalar of the intrinsic metric j ab of the 
brane world sheet is given by the Gauss formula 



R = K z ~ KlK h a = K z - Q , 



(32) 



where we introduced the notation Q for K^K^. Now we 
can rewrite S of (T25l) as 



with 



S = aK 2 



14 



bQ 



24 



24 



(33) 



(34) 
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and the curvature scalars K and Q are 



K 



1 

F 
1 

F 2 



r_0£ 
2 



B 
2F 2 



r 2 e 2 j 2 



B 2 
4i?4 



(n + l)f6 
fOB 



n cot 9 



(35) 



+/ 2 +' 



F 2 



cote 



where 



\/l + fr 2 6 , 
B = rf + 2/1 (9 + 4r/fl . 



(36) 



(37) 
(38) 



With the obtained expressions for the curvature scalars 
we will write up the Euler-Lagrange equation of the cur- 
vature corrected problem, and look for its regular solu- 
tions within a pcrturbative approach. 



IV. THE BRANE EQUATION 

Since the thickness corrected effective action ([2TJ)) also 
depends on the second derivative of 9 (as one can expect 
it from the presence of the curvature terms), the dynam- 
ics of the perturbed brane is described by the 4i/i-order 
Euler-Lagrange equation, 





i 9C ) 


d 


fdC\ 


8C 


dr 2 




dr 


\d'e) 


1 do 



= 



(39) 



Plugging the effective Lagrangian (f24|) into (|39|) we get 



d_ Cp 
d^ {jjf 

( 8(£ 6) 
dr 2 V 89 



8Cp 
89 



d 
dr 



8(£ 5)\ t 8(C S) 
86 



00 



.(40) 



At this point, to solve (|40|) . we choose a different path 
than the one that was followed by Frolov and Gorbonos 
in [2l[. First, we do not restrict our investigations to the 
Rindlcr zone and/or to the near critical solution of the 
thin problem; second, we follow more closely the treat- 
ment of Carter and Gregory in [20| . in the sense that, 
since the curvature corrections in (|23[) are second order in 
the perturbation parameter £/L, we are looking for a per- 
turbed solution whose perturbation is also quadratic in 
the same parameter. Hence we shall proceed by express- 
ing the solution of (|40|) as a power series in e = £ 2 /L 2 
as 



(r) = 9(r) + eip{r) + o(e 2 ), 



(41) 



where 9 is the solution of the zero thickness case, ip de- 
scribes its thickness correction at order of e, and we ne- 
glect all the higher order terms. 



Since the curvature scalars get larger and larger as we 
get closer to the singularity, the corresponding dynamical 
length scales become smaller and smaller. In the case of a 
fixed brane thickness £, this implies that the perturbation 
parameter e becomes larger too. To ensure the validity of 
the present perturbation method, we have to guarantee 
the condition that eip(r) <C 9(r), with the magnitude of 



e 4 < o.oi 



with respect to the second order pcrturbative treatment. 
As a consequence, for every fixed brane thickness value t, 
there is a minimum length scale parameter L that defines 
a condition on how close we can get to the singularity 
before the validity of the present perturbation formalism 
breaks down. 

By reducing the brane thickness £, of course, one can 
approach the singularity arbitrarily close. In this work 
however, we will chose a different approach. We fix the 
value of the perturbation parameter e = 10~ 2 . This 
means that we consider the thickest brane configurations 
that are allowed by the present perturbation method. 
This choice imposes the condition 



% <1 



and since < 9 < ^ , the corresponding condition on the 
amplitude of the perturbation function is 

7T 

M<2 ■ 

The above requirement tells us how close we can go to 
the singularity (i.e. how small the dynamical length scale 
L can be) before the present perturbation method breaks 
down for the case of the thickest possible brane configu- 
rations. 

Let us now substitute ([4*T|) into ([4*0]) , and after perform- 
ing the series expansion up to s, we obtain the following 
linear, inhomogencous equation for the perturbation 



where 



ip + qitp + q (p + q = 0, 



39 2 B 3 



qo = 



sin 



29B 2 
1 

A 9 {4) + A!(9 (3) 



A 2 9 +A 3 9 
A 8 f 2 + A 9 f + A 



10 



(42) 

(43) 
(44) 

A 4 / (3) 
(45) 



1 ~ 8r 4 /_F« 

+ A 5 f + A 6 f + A 7 f 



The coefficient functions Bi are the same of (fT3|) . while 
the coefficient functions Ai are given in the Appendix. 

Hence, the perturbed problem reduces to solving a 
second-order, linear equation for ip with a very compli- 
cated source term. In the next section we examine the 
asymptotic behavior and the regularity conditions of (|42p 
near the horizon, and in Sec. VII we discuss its solution. 
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V. REGULARITY CONDITIONS 

For a brane crossing the horizon (black hole embedding 
branch), the unperturbed solution 9 has the expansion 
(fH|) near ?' . hi this limit the source q has an a(r)/ f sin- 
gular behavior as r — > r$ with a regular c«o = ct(ro). The 
explicit expression for a is given in the Appendix. Hence, 
in the black hole embedding (supercritical) case the per- 
turbation equation (|4"2"|) has a regular singular point on 
the horizon. Requiring the regularity for tp at the horizon 
implies the boundary condition 



1 

m 



sin 2 6q 



a 



(46) 



where m = N — 3. Thus the perturbations for the super- 
critical solutions are uniquely determined by the initial 
value <p . For the later purpose of making a direct com- 
parison with the thin solution, we shall require the same 
boundary conditions for the perturbed and unperturbed 
brane solutions, therefore we choose tp to be zero and 
hence we get 



up 
in 



(47) 



If the brane does not cross the horizon (Minkowski 
embedding branch) the brane surface reaches its minimal 
distance to the black hole at r = r\ > r$. Near this point 
9 has the asymptotic form (|15|) . and one can find that q 
has a singular behavior 



Cl 



C3 



C5 



as r 



C5 



y/r — r\ (r — ri) 3 / 2 (r 
ri, where 



n) 5 /2' 



n(n— 
8(n+l)ri 



2(a + 2b){n + 1)/ + (a + 3b)rf 



(48) 



(49) 



and the explicit forms of c\ and C3 are given in the Ap- 
pendix. Thus the perturbation equation (|42|) . in the sub- 
critical case, near the vertical axis (9 = 0) has the asymp- 
totic form 



<P - 



n + 3 . 



2(r-n) 

Cl 



C3 



4(r — ri) 2 r — n 

C5 



1/r — r[ (r — ri) 3 / 2 (r — ri) 5 / 2 



= 0, 



where 



nr\ no 



(50) 



(51) 



rf?7 2 /(n) 12 77 

Equation ([50)) shows that no regular solution of 
exists at finite r\, unless C5 disappears. From (|4"9")) we can 
see that there are two possible ways to satisfy this condi- 
tion. One way is to choose the coefficients a and b in the 
appropriate way to make C5 disappear. Even though in 
([51)) we have already fixed the values of these parameters 



(we identified them with the constants that are obtained 
in the relativistic domain wall case in [20l]). however, so 
far we have not used their explicit values. Thus we could 
still consider the general case with arbitrary a and b. The 
problem with this way of regularization is that the con- 
stants a and b would be strongly dependent on the initial 
condition parameter r\, and hence provide a highly un- 
stable solution. 

The other way of making C5 disappear is to choose 
the brane dimension parameter n ■ — or n = 1. In- 
deed, instead of the previous option, we will follow this 
direction in the regularization process of (|50[) and keep 
the parameters a and b fixed as given in (|34p . This is 
also in agreement with what was found numerically in 
[2lj . (It can also be expected that, for n > 1, when 
additional symmetries in the angular directions are im- 
posed, the n = 1 case should be sufficient to discuss the 
brane configurations. Obviously, this is not the most 
general solution, and a different approach should be used 
to discuss the general case). The reason of the irregular 
behavior of the perturbations for n > 1 is due to the 
asymmetry between the sub- and supercritical configura- 
tions induced by the curvature corrections. It also seems 
obvious that going beyond this approximation, within a 
perturbative approach, will not resolve the problem. In 
[211 1 the authors suggest that quantum corrections may 
cure the above pathological behavior. Here, we take the 
simpler viewpoint by recording the fact that since the 
thin solution is not smooth on the axis, the perturbation 
method we used breaks down near this region for con- 
figurations outside the horizon. This, however, does not 
mean that physically reasonable solutions do not exist, 
as the explicitly constructed field theoretical domain wall 
solutions clearly show [25|, [26|. Hence, in the geometrical 
Dirac-Nambu-Goto approach, the thick solutions appear 
to deviate significantly from the thin ones near the axis of 
the system, and thus, our perturbative approach around 
the thin solution can not provide regular thick solutions 
for n > 1. In order to study the subcritical solutions 
and the thickness corrections to the phase transition in 
the general case, one needs to solve equation (|40[) in a 
non-perturbative way. In the following section we write 
down the explicit Euler-Lagrange equation for the non- 
perturbative problem, but due to its very complicated 
form we do not discuss its solution in the present pa- 
per. We rather stay with our perturbative approach and 
complete its full analysis in the next coming sections. 

In the following we provide the results for those cases 
where our perturbative approach leads to completely reg- 
ular solutions: 1-dimcnsional string or a 2-dimcnsional 
sheet, where n = D — 2, and D is the total number of di- 
mensions of the brane spacetime. In both cases, from the 
requirement of regularity, we get the asymptotic behavior 
of f(r) near the point n as 



ip = «V r — 7*1 + p(r - ri) 3 / 2 + 



(52) 
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where 



VII. THE PERTURBED SOLUTION 



2c 3 



(53) 



Thus the perturbations for the subcritical solutions are 
uniquely determined by the value of the parameter r\ . 



VI. THE GENERAL EULER-LAGRANGE 
EQUATION IN THE SUBCRITICAL CASE 

In the previous section we discussed that the pertur- 
bation method (|4Tj) breaks down around 9 = (the ver- 
tical axis) for cases n > 1 in the Minkowski embeddings, 
and hence one has to consider the full 4t/i-order Eulcr- 
Lagrangc equation (|40|) in searching for a regular solution 
of the curvature corrected problem. The equation of mo- 
tion in this general case is a very complicated and highly 
nonlinear one, 



where 



2i = 



Ti(6, 9, 9, f, /, r)9^ 
T 2 (e,9,e,f^,f,f,f,r) 



irf + /(2(2 + n) + r9(2ncot ( 
nr cot 99) — rf) 



(54) 



rfF 2 
+ r(0(2/(n-3 



(55) 



and T 2 is given in the Appendix. 

A detailed analysis of this equation, and a systematic, 
non-perturbative study to find its possible regular solu- 
tions would take us far from our essentially pcrturba- 
tive approach, and hence it lies beyond the scope of the 
present paper. Nevertheless, it does not seem completely 
impossible to perform this analysis, and to make a first 
step in this direction, we write up the asymptotic behav- 
ior of (|54|) near 9 = 0. This equation (obtained by taking 
the series expansion of (f5l)) around 9 = up to linear 
order) reads 



S3 

n3 



q2 



Si 



S + S9 = , 



(56) 



where the functions S, So, Si, S2, and S3 are given in 
the Appendix. With taking a look on the explicit forms 
of the Si functions, one can see that equation (|56p is 
essentially just as complicated as ([M]). and hence we will 
not discuss its properties further here. It is important to 
note however, that the study of this problem and finding 
its regular solutions are crucial for the understanding of 
the complete picture behind the phase transition in the 
curvature corrected system. 



A. Far distance solution 

As r — > 00 we can use the asymptotic form of 9 given in 
(fTS]) and (| 19|) . In the case of supercritical embedding, if 
n > 1 it is enough to consider the leading order p/r term 
to obtain the asymptotic behavior of the perturbation 
equation (14"2")) and we get 



n E 

^ + ^5 =° - 



where 



E 



-4p [(a + 2b) (n - 2) 2 - a(n - 4)] 



(57) 



(58) 



Equation (|57p can be integrated in a closed form, and its 
solution reads 



if - P + P ' 



2(n-3)r 3 



(59) 



From l|59p we can see that in the case of n — 3 the 
equation (|57p develops a resonance source term and hence 
(|59]l is not a solution. In this special case the solution 
reads 



P P' 
^=7 + 73 



£[l + 21nr] 

4r 3 



(60) 



In the case of n = 1 the perturbation equation (|4"2"|) , at 
far distances takes the asymptotic form 



1 
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-ip- 



-ip- 



E t + E 2 In r 







where 



E 1 = 4[a(llp'-4p) + 6(Tp'-2p)] J 
E 2 = -8p'(2a + b). 



E 1 +E 2 (l+hir) 
4 r 3 



The solution of (JeTJ) is 

P + P' lnr 

<P = 



(61) 



(62) 
(63) 



(64) 



and the integration constants P and P' in (|59|) . (|60| and 
(|64p are well defined continuous functions of tp or r\. 



B. Perturbations near the horizon 

After the analysis of the regularity conditions and the 
asymptotic behavior of the perturbation equation near 
the event horizon, it is now fairly straightforward to solve 
P2"| numerically. There is, however, one last thing that 
we have to discuss before presenting the numerical re- 
sults. In Sec. IV we have already mentioned that the dy- 
namical length scale of the problem has to be identified, 
in order to know how close one can get to the singularity 
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before the present pcrturbativc approach breaks down. 
Since 

y ~ max{A, y/\R\}, (65) 

we can calculate L for every initial parameter 9 Q and 
r\ from the supercritical and subcritical thin solutions 
respectively. We can then plug L into the curvature 
coefficients a and b of (|33| - ((34|) . and solve the pertur- 
bation equation. The requirement of how small 9o can 
be, i.e. how close we can get to the singularity before 
the quadratic perturbation approximation breaks down, 
is then obtained from the condition that \<p\ < This 
condition uniquely determines the length scale of the per- 
turbation problem in every chosen bulk (N) and brane 
(D) dimensions. 

In TABLE I we have listed the approximate minimal 
values of the initial parameters 8q of those thick brane 
configurations that can still be addressed by the present 
perturbation method. These values are used in the nu- 
merical calculations to obtain the dynamical length scale 
of the perturbed system in various bulk and brane di- 
mensions. It is interesting to notice that with increasing 
co-dimension (N—D) the minimum value of 9q is decreas- 
ing. Another tendency is that keeping the co-dimension 
fixed while increasing the bulk dimension, the minimum 
8q is also increasing. 



D\N 


5 


6 


7 


8 


9 


3 


6.82° 


4.52° 


3.4° 


2.7° 


2.25° 


4 


10.29° 


7.14° 


5.43° 


4.34° 


3.64° 


5 




9.85° 


7.56° 


6.18° 


5.14° 


6 






10.12° 


8.26° 


6.96° 


7 








10.48° 


8.82° 



TABLE I: The approximate values of the initial parameter 
6o for N — 5, . . . , 9 and D = 3, . . . , 7, where the present 
perturbation method reaches its limitation, i.e. where \tp\ ~ J. 



In the rest of this section we present the numerical solu- 
tion of the perturbation equation (|42[) in various dimen- 
sions. First we consider the case when n > 1, i.e. when 
there is no regular subcritical perturbations. We analyze 
this type of solution in the case of n = 2 with varying 
bulk dimensions. Later on we consider the n = 1 case, 
both for the subcritical and supercritical solutions. 

1. The n = 2, supercritical case. 

In FIG. 2 we have plotted the perturbations of a D = 4 
dimensional brane embedded into a N = 5 dimensional 
bulk spacetime. By varying the initial parameter 9q that 
describes the brane position on the horizon [6 measures 
the inclination from the vertical axis) until we get to its 




FIG. 2: The picture shows the numerical solution of the per- 
turbation equation (|42p with varying the parameter 8o of the 
brane initial inclination in the case of iV = 5, n = 2, black 
hole embedding. 



minimum, one can find four regions with qualitatively dif- 
ferent behavior. In the first region (top left), the pertur- 
bations are positive from the horizon to infinity, and their 
maximum amplitude grows moderately as 9q decreases. 
In the second region (top right), the perturbations are 
negative in the vicinity of the horizon, but soon become 
positive and after some increasing phase they start de- 
caying. In the 3rd region (bottom left), the picture is 
similar to the 1st one with the difference that the ampli- 
tude starts growing more rapidly as 9q decreases. Finally, 
in the Ath region (bottom right), the perturbations are 
positive in the vicinity of the horizon, and there is a sign 
change near the horizon before they start decaying. 

In FIG. 3 we have plotted the perturbations of the 
near minimal 9q zone (\ip\ < m different bulk dimen- 
sions. By increasing the bulk dimension and keeping the 
brane dimension fixed, we found that the perturbations 
amplitude are decreasing and hence the minimal 8q is 
also decreasing. This means that we can get closer to the 
singularity 9 = before the perturbation approximation 
breaks down, i.e. \ip\ ~ |. It is also interesting to notice 
that the sign change, in the near singularity region so- 
lutions, gets closer and closer to the event horizon with 
increasing bulk dimensions. 

In FIG. 4 we have plotted the corresponding thick 
(perturbed) and thin (unperturbed) brane configurations 
with the initial parameters 9q = \ (bottom curves) and 
y^j- (top curves) respectively, in a cylindrical coordinate 
system with N = 5 bulk dimensions. On the picture, the 
black curve represents the black hole's event horizon, the 
blue thin curves the unperturbed brane solutions, while 
the red, thick curves the perturbed ones. 
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FIG. 3: The picture shows the numerical solutions of the 
perturbation equation (|42[) in the case of a D = 4 dimensional 
(n = 2) brane embedded into a iV = 6, 7, 8, 9 (top left, top 
right, bottom left, bottom right, respectively) dimensional 
bulk spacetime. The graphs belong to the initial conditions 
that are close to the minimal 9o, i.e. when the perturbation 
method reaches its limit \(p\ ~ -|. 




FIG. 4: The picture shows the thick (red) brane configura- 
tions together with their thin (blue) counterparts in a cylin- 
drical coordinate system, in the case of an N — 5, n = 2 black 
hole embedding. The initial conditions are #o — f (bottom 
curves) and ^f-g (top curves), and the thickness parameter £ is 
chosen to be large for the purpose of making the effects visible. 
The black curve represents the black hole's event horizon. 



2. The n — 1 case. 

In the case of n = 1, i.e. a D = 3 dimensional brane, 
the perturbation equation (|42|) has regular solutions both 
in the supercritical and the subcritical cases. For the 
black hole embedding branch the perturbations are very 
similar to the n = 2 case. The only remarkable differ- 
ence, compared to the results that are plotted in FIG. 2, 
is that in the n = 1 case the 2nd region of FIG. 2 is 
missing, i.e. the perturbations are always positive in the 
very vicinity of the horizon. 



The numerical solutions of the perturbation equation 
in the case of the Minkowski embedding branch can be 
seen in FIG. 5. We have chosen two regions, charac- 
terized by their minimal distances from the black hole's 
event horizon. The first picture (top left) in FIG. 5 shows 
the perturbations in cases when the brane is very close 
to the horizon, 1.001 < r\ < 1.01. The second picture 
(top right), however, belongs to the solutions that are 
a bit further away from the horizon, in the region of 
1-1 < T\ < 3. The top pictures belong to N = 5, whilst 
the bottom pictures to N = 6 embeddings. 



<fi if 




FIG. 5: The hgure shows the numerical solutions of the per- 
turbation equation (|42|) in the cases when a D — 3 dimen- 
sional brane is embedded into an N = 5 (top pictures) and 
an N = 6 (bottom pictures) dimensional bulk. The pic- 
tures belong to the Minkowski embedding branch with initial 
parameter regions 1.001 < n < 1.01 (left hand side), and 
1.1 < ri < 3, (right hand side). 

From FIG. 5 we can see that the perturbations in the 
near horizon zone start with a negative maximum and 
soon after there is a sign change before they begin to 
decay. The far zone solutions start also with a nega- 
tive maximum but they begin to decay without changing 
their sign. From the top to the bottom, we have plotted 
the same solutions changing the bulk dimension from 5 
to 6. In accordance with the corresponding black hole 
embedding solutions, we can see that the perturbations' 
amplitudes are also decreasing with increasing bulk di- 
mensions. 

In FIG. 6 we have plotted two corresponding thick 
(perturbed) and thin (unperturbed) brane configurations 
with initial parameters r\ = 1.001 (bottom curves on the 
Z axis) and r\ = 1.04 (top curves on the Z axis) re- 
spectively, in a cylindrical coordinate system with N = 5 
bulk dimensions. On the picture, the black curve rep- 
resents the black hole's event horizon, the blue curves 
the unperturbed- (thin), while the red curves the per- 
turbed (thick) brane solutions. We can see on the picture 
that the near horizon perturbed configuration intersects 
its thin counterpart as the perturbation changes its sign 
before it starts decaying. There is no intersection, how- 
ever between the thick and thin configurations at the top 
curves, since the perturbations in the far zone do not 
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FIG. 6: The picture shows the thick (red) brane configura- 
tions together with their thin (blue) counterparts in a cylin- 
drical coordinate system, in the case of an N = 5, n = 1 
Minkowski embedding. The initial parameters are n = 1.001 
(bottom curves) and rj = 1.04 (top curves), and the thickness 
parameter £ is chosen to be large for the purpose of making 
the effects visual. The black thick curve represents the black 
hole's event horizon. 

change their sign before they start decaying. 

VIII. ENERGY PROPERTIES AND PHASE 
TRANSITION 

In obtaining the solutions of the curvature corrected 
brane action with a second order perturbative ap- 
proach in the thickness, we concluded that our pertur- 
bative method breaks down in the vicinity of the axis 
of the system, for the subcritical solutions with n > 1, 
and hence one needs to find a new, exact solution for a 
general discussion. Nevertheless, we also found that in 
the special cases of n — and n = 1, the perturbation 
equation (|42|) can be regularly solved in the subcritical 
case too, whereas we also showed that the supercritical 
solutions are regular and can be solved without any dif- 
ficulty for any n. In the following we will discuss the 
phase transition between the Minkowski embedding and 
black hole embedding topology of the thickness corrected 
brane solutions with n = 0, 1. Since we are interested in 
brane solutions we present the results of the n = 1 case. 

For the discussion of the properties of this possible 
transition we use the method presented by Flachi et. al 
in (22| , where a quasi-static evolution of a brane has been 
considered from the equatorial configuration. In their 
work, the authors demonstrated that if one plots the dif- 
ference in the energy AE between a brane configuration 
(that is quasi-statically evolved from the equatorial con- 
figuration and uniquely defined with its initial parameter 
g = 9q or n) and the equatorial brane configuration, with 
respect to the cylindrical distance parameter 

Zoo(g) =Too cos 6(g, r^, 

then, the obtained plot exhibits a loop, i.e. an instability 
zone, that is usually a typical sign of a first order phase 
transition in dynamical systems. The parameter is a 
practically chosen maximum radial distance value, until 
the equations are integrated numerically from the hori- 
zon. It is so large that it has the property that the near 



AE AE 




FIG. 7: The figure shows the energy difference of a brane con- 
figuration (that quasi-statically evolves from the equatorial 
configuration) with respect to the equatorial configuration, 
as a function of the parameter in the case of N = 5 bulk 
dimensions. The top pictures belong to the thin, while the 
bottom pictures to the thick system. The red (blue) curves 
represent the black hole (Minkowskian) embedding branch. 
The pictures on the right hand side are the zooms into the 
near region of Omin , where the effects of the perturbations are 
the largest, and the difference become apparent between the 
thin (top) and thick (bottom) cases. 




1 2 3 4 5 6 Z ™ 4 ^.60 3.65 3.70 3.75 3.8?' 



AE AE 




FIG. 8: The pictures on this figure are the corresponding ones 
to FIG. 7, in the case of TV = 6 bulk dimensions. 

horizon behavior of the function AE(Z 00 ) is essentially 
not affected by any further increase in r^. 

The energy of a brane with a given initial parameter g 
can be calculated as 

E(g) = - [°° L(g,r) dr , 

and 

AE(g) = E(g) - 25(f) . 

In FIG. 7 and FIG. 8 we have plotted the curves of 
AE(g)[Z OQ (g)] with varying g = 6q in the region of 
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(fi^mm) for the black hole embedding branch (red 
curves), and with varying g = r\ in the region of 
(1.001,10) for the Minkowski embedding branch (blue 
curves), in N = 5 and 6 bulk dimensions respectively. In 
both figures, the top left picture represents the unper- 
turbed, thin brane solutions, and the bottom left picture 
the thick ones. In the plots it is easy to see the pres- 
ence of the loops which are the sign of a first order phase 
transition. This implies that small thickness perturba- 
tions induced by curvature corrections in the effective 
brane action do not change the qualitative behavior of 
the transition in the dimension where a regular pertur- 
bative solution exists. 

To make the effects of the corrections visible, we en- 
larged the region where the perturbations are getting 
larger, i.e. near the singularity. We zoomed into the close 
neighborhood of the initial parameter 9 m i n (listed in TA- 
BLE. I) for both the thin and thick curves. The results 
can be seen on the right hand side pictures of both figures. 
As a conclusion we can say that within the perturbative 
approach, small thickness perturbations induced by cur- 
vature corrections in the effective brane action have small 
effects on the AE(Z 00 ) plots of the system, and they do 
not change the qualitative behavior of the phase transi- 
tion. 

It is very important to note however, that the obtained 
results are valid only for the exceptional cases where the 
problem can be regularly solved within the perturbative 
approach. In principle, one can expect a different solu- 
tion to the thick brane equation with a non-perturbativc 
treatment in the subcritical cases, thus the behavior of 
the phase transition can also be expected to be different 
from the picture above. 



IX. CONCLUSIONS 

In this paper we studied thickness perturbations to 
static, D-dimensional, Dirac-Nambu-Goto branes em- 
bedded into higher dimensional, spherically symmetric, 
black hole spacetimes. The perturbations originate from 
higher order, curvature corrections added to the thin 
brane action (l9l . [2p| . and are quadratic in the thick- 
ness of the brane. Following the treatment of [2(|, and 
applying a linear perturbation method with the pertur- 
bation parameter e = £ 2 / L 2 , we derived the general form 
of the perturbation equation for a brane that is axisym- 
metric and has a form of a (D — l)-dimcnsional plane at 
asymptotic infinity. 

From the analysis of the asymptotic behavior of the 
perturbation equation, we found that there is no regu- 
lar solution of the perturbed problem in the Minkowski 
embedding case, unless the brane is a string, or a sheet. 
This restriction, however does not hold for the black hole 
embedding solutions, which are always regular within our 
perturbative approach. 

From the above results and also from the similar find- 
ings of [2l[ , we concluded that the absence of regular solu- 



tions above the dimension D — 3 implies, that the prob- 
lem can not be solved within perturbative approaches 
around the thin solution which is not smooth on the axis 
of the system. Hence, for a general discussion, one needs 
to find a new, exact solution of the curvature corrected 
problem, that is expected to behave differently from the 
thin solution with being smooth on the axis. Finding 
this new solution however seems to be a difficult task, as 
the equation of motion is a very complicated and highly 
nonlinear one. 

After the above conclusions, we provided the solution 
of the perturbation equation for various brane (D) and 
bulk (N) dimensions. The far distance equations are 
integrated analytically, while the near horizon solutions 
are obtained by numerical computations. The deforma- 
tions of the perturbed brane configurations are plotted 
and a comparison is made with the corresponding thin 
brane configurations with identical boundary conditions, 
for both types of solution. 

One motivation of this paper was to consider the effects 
of higher derivative, curvature corrections on the first or- 
der phase transition between the Minkowski and black 
hole branch, that is present in the unperturbed system 
El GJ, H3- With the solution of the perturbed prob- 
lem we found that within a perturbative approach, one 
can consider a phase transition between the two branches 
only in the cases of D = 2 or 3. We investigated the prop- 
erties of this transition in the case of a D = 3 brane, and 
found that small thickness perturbations do not mod- 
ify the qualitative behavior of the phase transition, i.e. it 
remains a first order one, just as in the case of zero thick- 
ness. 

Since our perturbative approach does not provide a 
regular thick brane solution for dimensions D > 3, we 
cannot answer, in the most general way, the question 
whether higher order, curvature corrections in the effec- 
tive brane action can change the order of the phase transi- 
tion in the BBH system or not. Although we expect that 
small corrections may not change the picture too much, 
as they are quadratic in the thickness of the brane, how- 
ever a definitive answer can only be given after a new, 
non-perturbativc solution has been found. Further study 
to address this question is in progress and we hope to 
report on it in a forthcoming publication. 
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A = -6A(a + b)r 4 f 2 F 4 , (A.l) 

A x = -6A(a + b)fFr 3 csc9{29 fFnr cos9 + (Afr + f (2(2 + n) + 9r 2 (9{2f(n - 3) - 10fr9 - fr)))) 81119), (A.2) 

A 2 = -16O(a + 6)/ 3 r 6 (60 2 /r 2 - 1), (A.3) 

A 3 = 48(a + 6)/ 2 r 4 esc 9(9r(69fnr cos + 80 3 / 2 w 3 cos 9 - 59 2 fr 2 (fr - 2/(n - 3)) sin 9 

+ 10(/(4 + n) + 3/r)sin0) - 2ncos0), (A.4) 

A± = -l69fF 2 r 4 (2a + b + a9 2 fr 2 ), (A.5) 

A 5 = AFr 2 esc 0(4an cos + r(-20 2 /nr(7a + 46 + 2 /r 2 (13a + 46 + 4a0 2 /r 2 )) cos 9 

+ (2a0 7 / 4 (l + n)r 6 - 8(2a + 6)/0r + 4 /V(0(8a - 66 - 6an) + 8a0r + (a + 6)/0 3 r 3 ) 

+ 40 2 /V(0(2a + b- 2(4a + 6)n) + (7a + 66)0r + (2a + 6)/0 3 r 3 ) + 4/(-20(3a(2 + n) + 6(5 + n)) 

- (lOa + 96)0V + (4a + 56)/0 3 r 3 ))sin0)), (A.6) 

A 6 = 26>V(80(-71(a + 6) + 45an + 426n + 2 (a - 26)n 2 )r + 40(6(-8O + (69 - 7n)n) 

+ a(-80 + n(54 + (29 - 3n)n)) + 20(a(37 - 17n) + 6(25 - 8n))nr 2 cot 9) - 40* w 3 cot 9(2n(-b + an) 

+ 3/(a + 26 + 3an)r + 2(-2 + n)(2a + 6 - an) esc 2 0) + /0 5 r 5 ((3a + 6)/(l + n)r 

+ 4an(n - (-1 + ra) esc 2 0)) + 40 r 2 (/(26(8 + n) + a(13 + ?i(27 + 5n)))r + rt(—a(— 25 + n)n 

+ 6(-3 + lln) + (6(13 - lln) + a(28 + (-27 + n)n)) esc 2 0)) + 40 2 r(?i(6(57 - 31n) 

+ a(89 - n(26 + 19n))) cot + 0r 2 (3/(a + 26 + 3an)r + 2n(-6 + 5an + 46n 

+ (6a + 56-5a7i-46n)csc 2 0)))), (A.7) 

A 7 = 2 r 3 (-80(a(-24O + n(-68 + 13n)) + 6(-240 + n(-47 + 16n)))r - (a + 6)./ 3 7 r 9 
+ 2(5a + b)f 2 9 6 nr 7 cot + 80(6(48 - 13(-4 + n)n) + a(48 + n(60 + (8 - 5n)n)) 

- 20n(166n + a(-9 + 34n))r 2 cot 0) + 80 2 r(n(6(42 - 43n) + a(78 - n(76 + 15n))) cot 

+ 0r 2 (16(a + 6)n + ./(-14a- 236 + 43an+ 166n)r + n(16a+ 196- 13an - 166n) cot 2 )) 

+ 20 4 r 3 (-(5a + 6)0./ 2 ?- 4 + 4ncot0(n(76- 5an) - /(7a + 146 + 42an + 46n)r 

+ (n - 2)(-12a - 76 + 5an) esc 2 0)) + 2/0 5 r 5 (./(56(2 + n) + a(14 + 25n))r + 2n(-b + Uan + Abn 

+ (12a + 56 - lion - 46n) esc 2 0)) + 40 3 ?' 2 (/r(a(76 + (150 - 19n)n) + 6(58 + (9 - An)n) 
+ 6(3a + 46)0w 2 cot 0) + 2n(a(44 - 15n)n + 6(-6 + 17n) + (6(20 - 17n) 

+ a(50 + n(-62+ 15?i)))csc 2 0))), (A.8) 

A s = -2r(80(a(2 + n)(4 + 5n) + 6(8 + n(15 + 4n)))r + 3(2a + 6)/ 3 7 r 9 - 23(a + b)j 2 9 & nr 7 cot 

+ 80n(n(36 + a(2 + n)) + 0(28a + 256 + 17an + 86n)r 2 cot 0) + 2/0 5 r 5 (-3/(56(-l + n) + a (4 + lln))r 

+ (a + 46)n(-2 + n + ncos[20])csc 2 0) + 20 4 r 3 (3(7a+ lO6)0/ 2 r 4 + 2ncot0(3n(-36 + an) 

+ /(74a + 476 + 53an + 8bn)r + 3(-2 + n)(4a + 36 - an) esc 2 0)) + 40 3 r 2 (/r(2(6(-78 + n(13 + 2n)) 

+ 2a(-42 + n(-5 + 6n))) + (47a + 266)0w 2 cot 0) + n(6(3 - 9n) + an(-21 + 19n) 

- (n - l)(-96 + a(-24 + 19n)) esc 2 0)) + 40 2 r(n(216n + a(-6 + n(50 + n))) cot 

+ 20r 2 (3/(a(-48 + n) + 6(-46 + 3n))r - (a - 26)n(-2 + n + n cos [20]) esc 2 0))), (A.9) 
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A 9 = 2(-8(6 - o(-2 + n))(-l + n)ncot6» + r(-9(2o + 36)/ 3 6> 5 r 7 + 2(-17o + 46)/ 2 6> 4 ?7r 5 cot (9 

+ A6(2fr(-2a(2 + ra)(2 + 3n) - 6(8 + n(15 + ra)) - (28a + 19b)6nr 2 cot 9) 

+ n 2 (-(b - 7an) cot 9 2 - 8a esc 2 9)) - 4flr(2/(20a(2 + n) + 6(40 + YIn))r + n((b - 7an) cot 6 2 + 8a esc 2 9)) 

+ 2<9 2 r(3(46a + 436)0/V + 2n cot 0(ri(56 + an) - /(74a + 416 + 4(4a + b)n)r 

- (n - 2) (56 + a(4 + ra)) esc 2 0)) + 2/0 3 r 3 (/(2a(62 + n) + 6(102 + 7n))r + n(3b - 2\an + Abn 

+ (6-46n + 3a(-8 + 7?i))csc 2 0)))), (A.10) 



A w = 4(2(10a + 76)/ 3 6» 3 r 6 + / 2 r 3 (-20(lOa(2 + n) + 36(6 + n)) - 2(18a + 176)0r - (26a + 76)0 2 ?ir cot 0) 

— n(6 + art) (—4 + n + n cos [20]) cot 9 esc 2 (9 + fr(9 9 f(l + n)(3(a + 6) + 5an)r 9 - 4(6 - 2a(-l + n))ncot 9 

- 29nr((b - 7an) cot 9 2 + 8a esc 2 9)) + 9 & f(a + 6 + an)r 7 (20(-13 + 7n)r + 26 (-24 + n(17 + n) 

+ 0r((13-7n)ncot0 + 0r(/(-l + n 2 ) +2n(n- (-l + n)csc 2 0)))))). (A.ll) 



a = - -^-[-2(lOa + 76)r 6 / 3 + 2n cot 6((b + an)((n - 2) esc 2 9 - n) - 2ar 2 /) + 2r/(2n(2a + 6 - 2an) cot 



1 

2T A 

+ r'e(n((b - 7an) cot 2 9 + 8a esc 2 9) + 4(2a + 6)r 2 /)) + r 3 ./ 2 (0(2Oa(2 + n) + 66(6 + n) + (26a + 76>w cot 99) 
+ 2(18a+ 176)r0)) ] (A.12) 



ci = (l/(960r; 5 / 2 7-i))(3840(6(180 + n(109 + 8n)) + o(180 + n(100 + 17n))) + r)n(15fr)n (64(6(272 

+ n(164 + 13n)) + a(272 + n(166 + n(34 + n)))) + 16/?/ 2 (a(39 + 5n(4 + n)) + 6(69 + n(39 + 4n)))r 2 

- 4/V(6(4 + n) + a(14 + bn))r\ + (a + b)f 3 jfrf) + 120n 6 / 4 (a + 6 + on)(-l + n 2 )r?(-r/ + 9ricr) 
+ 4r? 3 / 3 r 2 (-4(a + 6 + an)(-15r/ 2 (-68 + n(23 + 9n)) + 5n 4 (14 - 5n)nri + r/ 6 (-l + n)nr 2 

+ 1577(39 + n(-79 + 28n))ricr - A5n 3 n(l + 2n)r\o + 180(-1 + n)nr 2 l( r 2 ) 

+ 15?7 3 (1 + n)r 2 (/(4?7(-l + n)(a + 6 + an) + 9(3(a + 6) + han)ria) + r 1 (2af {3) rjr 1 

+ f(3bn + a(577(l + n) + 18ricr))))) + 5/(64(6(6?7(436 + 3n(55 + 2n)) + ?/ 3 n(62 + 19n)n 

+ 3(2652 + n(1397 + 91n))ri<r) + 0(677(436 + ra(190 + n(15 + 2n))) + ifn{8 + n)(10 + 3n)n 

+ 3(2652 + n(1304 + n(181 + 3n)))ricj)) - ?y 2 r 2 (48/77(6(60 + n(35 + An)) + a(36 + n(32 + 5n)))n 

+ 3(a + 6)/ 3 77 4 r 4 (2?7 + 9ri<j) + ,/Vr 2 (9(a + 6) ./Vrf + 4(-677(5a + 76 + 2(8a + 6)71) + (5a + b)n 3 n ri 

- 3(5a + 6)(-3 + 5n)ria)) - 8/(3/?7 3 (6(5 + n) + a(17 + 5n))rf + 2(a(-6r/(-58 + n(16 + 25n)) 
+ 2?7 3 ?i(12 + 5n)n + 3(564 + (132 - 5n)n)ricr) + 6(-6t?(-94 + n(5 + 4n)) + 2t7 3 ti(13 + An)n 

+ 3(564 + (173 - An)n)na))))) + 2nf 2 r 1 (m(b(-30r l 2 (-l57 + n(3 + 14n)) + 20? ? 4 ?i(13 + ri)n - 7?7 6 nr 2 

+ 1577(1884 + (115 - 7\n)n)r x a + 5n 3 n(66 + 137i)r 2 cr + 45(1392 + n(527 + 13n))r 2 cr 2 ) 

+ a(-30?7 2 (-157 + n + 36?i 2 + 7?i 3 ) + 20t7 4 ti(27 - 2(-2 + n)n)n + n 6 n(-10 + n(9 + 4n))r 2 

+ 1577(1884 + n(276 + (-191 + n)n))na - 75?7 3 n(-14 + (-2 + ?i)n)r 2 cr 

+ 45(1392 + 7i(522 - (-19 + n)n))rla 2 )) + 77 2 r 2 (15/V(l + n)rf(r)(7b + a(9 + 10n)) + 9(3a + b)na) 

- 10777-1 (3a/ (3) ?7(32 + n(15 + n))n - 3 (a + b)f 2 ifr 3 + 2/(3677(3 + n)(8 + n) + a{3n(-2 + n(A3 + n(18 + n))) 

- 77 3 ?i(9 + 2n) n + 3(-24 + n(33 + 7n))na))) + 2/(15?7 3 r 3 ((a + ^f^r/n + /(??(6(3 + n) + o(7 + 5n)) 
+ 9(a + 6)ricr)) + 2(10677(1877(7 + n - n 2 ) + ?7 3 n(7 + 2n)n - 3(-18 + n(23 + 7n))na) 

+ a(90?7 2 (ll + n(20 + n - 2n 2 )) + 20n 4 n(2 + n(6 + n))n - rf{-l + n)nr\ 

- 30?7(-9 + n(-19 + n(45 + 7n)))r x o + A5rfn{l + 2n)r\a - 180(-1 + n)nr\a 2 ))))))) \ r (A.13) 
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c 3 = (l/(192r7 3 /r 4 ))(-192(6(20 + n(47 + 9n)) + a(20 + n(40 + n(15 + n)))) 

+ ?yri(24r/ 5 / 3 (a + & + on)(-l + n 2 )r 2 - 3,/Vi (16(6(4 + n)(3 + 4n) + o(12 + n(12 + 5n))) 

- 4(5a + 6)/r/ 2 (l + n)r? + (a + &)/V»i) + 4r) 2 fn(3r) 3 (l + n)r?(/(3(o + 6) + San) + 2o/n) 

+ 4(a + & + an)(3r/(13 + (5 - 6n)n) + J] 3 n(l + 2n)n - 21(— 1 + n)nr lC r)) - 2/(16(&(3t?(36 - 5(-5 + n)n) 

+ 77 3 ?i(2 + 5n)n + 3(96 + 23n(7 + n))n<r) + a(-3r?(l + n)(-36 + 5n 2 ) + ifn{6 + n(10 + 3n))n 

+ 3(96 + n(150 + (37 - 3?i)n))r lC r)) + n 2 r\ (12a/r?(8 + n(7 + - 3(3a + 6)/V(l + n)r^ 

+ 2/(-3(a + b)/r/ 3 r? + 2(6677(2 + n)(3 + n) + 0(677(1 + n)(3 + n(8 + n)) - n 3 n(l + 2n)n 



21(-l + n)nntr)))))))l 



(A.14) 



480(a + 6)r 4 e/ 2 F 2 



n(-8e(l + 2n)(& + an) + 20r 2 + 206* e/ 4 (-13 + 7n)(a + 6 + an)r 6 



+ 10er(-4afr + /(8a + 46 - 8an + (26a + 76)/6> 2 r 3 )) + 4(? 2 / 3 r 4 (20(9 2 £(6(-ll + 6n) 

+ a(-19 + ra(8 + 3n))) - 240(a + 6)# 2 £r 2 + 10fer(#(-37a - 256 + 17an + 86n) 

+ 8 (a + 6)6> (3) r) + 0V(4e(-&(l + 2n) + 2a(-l + n + n 2 )) + 5r 2 + 5er(3/(a + 26 + 3an) 



er 



+ 4a/r))) + /V(2O0 e(6(2 + 19n) + a(-2 + n(44 + 3n))) + 240(a + 1 

- 60(3a + 4,b)6f6 3 er 4 ~ 5(5a + b)f 2 f) 6 er 6 + 4O0er(0(28a + 256 + 17an + 86n) 

+ 8(a + 6)0 {3) r) + 40V 2 (2e(-5(& + 2bn) + a(-8 + n(7 + 6n))) + 15r 2 

+ 5er(/(46(2 + n) + a(4 + 33n)) + (13a + 46)/r))) + 2/(20e(6 - o(-2 + n))(-l + n) 

+ 0r 2 (2O(28a + 196)0/W 2 + 0(8e(-2(a + 6) + (a - 46)?i) + 30r 2 + 5er(2(7a + 46) /V 

+ /(786 + 86n + 20a(7 + 2n) - (9a+116)/0 2 r 3 )))))) ] , 



(A.15) 



Si = 



32(a + b)r 4 ef 2 F 2 



n(-4en(6 + an) + 4r 2 + 40 e/ 4 (-13 + 7n)(a + 6 + an)r b 

2 



2er(-4a/r + /(8a + 46 - 8an + (26a + 76)/0~r 3 )) + 40 /V(40 e(b(-U + 6n) 
a(-19 + n(8 + 3n))) - 48(a + b)Q 2 er 2 + 29er(6(-37a - 256 + 17cm + 86n) 
8 (a + 6)0 (3) r) + tr 2 {2en(-b + an) + r 2 + er(3/(a + 26 + 3an) + 4a/r))) 
f 2 r 2 (4,9 2 e(b(2 + 19n) + a(-2 + n(44 + 3n))) + 48(a + b)Q 2 er 2 
12(3a + 4b)6f8 3 er 4 - (5a + 6)/ 2 6 £r 6 + 80er(0(28a + 256 + 17an + 8bn) 
8(a + 6)0 (3) r) + 40 4 r 2 (en(-56 + San) + 3r 2 + er(/(46(2 + n) + a (4 + 33n)) 
(13a + 46)/r))) + 2/(4e(6 - a(-2 + n))(-l + n) + 6>r 2 (4(28a + 196)0/er 2 
6(-8bsn + 6r 2 + er(2(7a + 46) fr + /(786 + 86n + 20a(7 + 2n) 
(9a+116)/0 2 r 3 ))))))] , 



(A.16) 



1 



16(a + 6)r 3 / 2 



n(2/0(6 + a(8 - 7n))r + 40 f(n - l)(a + 6 + an)r 4 



+ /(20n(8a + 6 - 7an) + 20(8a + 6 - 7an)r + fO (7an + Abn - 8a - 56)r 3 ) 
+ 20 2 /V(20(6(4n - 5) + a(5n - 6))r + 6»(6(7n - 9) + a((27 - 5n)n - 24) 
+ a/0 2 (7i-l)r 3 ))) 



5 3 



1 



(a + 6)r 4 / 2 



n(n - 2)(1 + r 2 f8 2 )(b + an + 2(6 - o(n - 2))r 2 ftf) 



(A.17) 
(A.18) 
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T2 = w^h^ [ -^ ef{a+b+an){ - 1+n2y 

- 2e 2 /V(240(a + b)6 2 ee(-3 + n)r 2 - 480(a + b)f ''er 3 + 860 2 er(-7l{a + b) 

+ 45an + 42bn + 2(a - 26)n 2 + 24(a + b)6nr 2 cot 8) + 40 J e(&(-8O + (69 - 7n)n) 
+ a(-80 + n(54 + (29 - 3n)n)) + 2#(a(37 - 17n) + 6(25 - 8n))nr 2 cot 9) 
+ A8 6 nr 3 cot 8{2en(b - an) - r 2 + er(-3/(a + 2b + 3an) - Aafr) 

- 2e(6 — a(— 2 + n))(— 2 + n) csc 6> 2 ) + 48 r 2 (en(-a(-25 + n)n + 6(-3 + lln)) 

+ (5 + 4n)r 2 + er(/(26(8 + n) + a(13 + n(27 + 5n))) + r(/(5a - 36 - 2an) - 2a/ (3) r)) 

+ en(6(13 - lln) + a(28 + (-27 + n)n)) esc 8 2 ) + A8 4 r(en(b{57 - 31n) 

+ a(89 - n(26 + 19n))) cot 8 + 8r 2 (2en(5an + b(-l + An)) + r 2 + er(3/(a + 2b + San) 

+ Aafr) + 2e(a(6 - 5n) + 6(5 - 4n))n esc <9 2 )) + /#V((3a + 6)/e(l + n)r + 2(2aen 2 

+ r 2 + (a + b)fer 2 - 2ae(-l + n)n esc 8 2 ))) + 2/ 2 r(3(2a + b)f 3 8 7 £r 9 

- 23(a + b)f 2 8 6 enr 7 cot + 88er(a{2 + n)(4 + 5n) + 6(8 + n(15 + 4n)) 

+ 6{a + b)8nr 2 cot 8) + 88e{n 2 (3b + a(2 + n)) - 90(a + b)8 2 fr 5 + 8n{b{25 + 8n) 

+ a(28 + 17n))r 2 cot 8) - 2/eV(2(a + 4b)en 2 + 9r 2 + 3er(/(56(-l + n) + a(4 + lln)) 

+ 8(a + 6)/r) - 2(a + 46)e(-l + n)ncsc 8 2 ) + 48^ r 2 [en{b(3 - 9n) + an{-21 + 19n)) 

+ (-7 - 4n)r 2 + er(2r(/(6(9 + An) + 2a(5 + 7n)) + (5a + 26)/ (3 V) + /(2(6(-78 

+ n(13 + 2n)) + 2a(-42 + n(-5 + 6n))) + (47a + 26b)8nr 2 cot 8)) - e(-l + n)n(-96 

+ a(-24 + 19n)) esc 8 2 ) + 2^V(3(7a + 10b)8f 2 er 4 + 2enr(f(b(47 + 8n) + a(74 + 53n)) 

+ 4(5a + 26) fr) cot 9 + 6ncot 6>(en(-36 + an) + 2r 2 + e(36 - a(-4 + n))(-2 + n) csc6» 2 )) 

+ A8 2 r(en(21bn + a(-6 + n(50 + n))) cot 8 + 8r 2 (A(a - 26)en 2 - 3r 2 + 6er(/(a(-48 + n) 

+ 6(-46 + 3n)) + (a + b)fr) - 4(a - 26)e(-l + n)ncsc8 2 ))) - 4$ fr 7 (28e(a + b 

+ an)(-2A + n(17 + n)) + 2#e(-13 + 7n)(a + b + an)r - 2^en(-13 + 7n)(a + 6 

+ an)r cot 6> + ^ 3 r 2 (e(l + n)r(/(3(a + 6) + 5an) + 2a fr) + 2(2en(a + b + an) 

+ (1 + n)r 2 ~ 2e(-l + n)n(a + b + an) esc 8 2 ))) - 2/(-9(2a + 3b)f 3 8 5 er s 

- 8e(6 - a(-2 + n))(-l + n)n cot 6> + 2(-17a + Ab)f 2 8 4 £nr 6 cot (9 - 40r 2 (-r 2 

+ 2er-(/(20a(2 + n) + 6(40 + 17n)) + (10a + 96) fr) + en(b - 7an) cot 8 2 + 8aen esc 8 2 ) 

- A8r(sn 2 (-b + 7an) + (-2 - n)r 2 + 2er(r(2/(3a(2 + n) + 6(5 + n)) + (2a + 6)/ (3) r) 
+ /(2a(2 + n)(2 + 3n) + 6(8 + n(15 + n)) + (28a + 196)Snr 2 cot 8)) + e(b 

+ a(8 - 7n))?i 2 csc6i 2 ) + 2f8 £ '' r 4 (en(-21an + 6(3 + 4n)) + 7r 2 + er(/(2a(62 + n) 

+ 6(102 + 7n)) + 126/r) + en(6 - 46n + 3a(-8 + 7n)) cscd 2 ) - 28 2 r 2 (-3{A6a + A3b)8f 2 er A 

+ 2n cot 6»(-en(56 + an) + Ar 2 + er(/(74a + 416 + 4(4a + 6)n) + (5a + 46) fr) 

+ e(-2 + n)(56 + a(4 + n)) cscd 2 ))) + f 3 r 3 (-A80(a + b)8 2 8s(A + n)r 2 

- 160(a + b)8 A er 3 + (a + b)f 3 8 9 £ r 9 - 2(5a + b)f 2 8 % enr cot 8 

- 8e<9 2 CT-(6(240 + (47 - 16n)n) + a(240 + (68 - 13n)n) + 36(a + b)8nr 2 cot 8) 

+ 88 3 e(b(-A8 + 13(-4 + n)n) + a(2 + n)(-24 + n(-18 + 5n)) + 30(a + b)f fr 5 

+ 28n(16bn + a(-9 + 34n))r 2 cot 8) + 88* r{sn(b(- 42 + 43n) + a(-78 + n(76 + 15n))) cot 8 

- 8r 2 (16(a + b)en + 3r 2 + er(/(-14a - 236 + 43an + 166n) + 6(3a + 26)/r) 

+ s(6(19-16n) +a(16- 13n))ncot6» 2 )) + 2e 6 r 3 ((5a + 6)e/ 2 £r 4 + 4ncot6»(-76en + 5aen 2 

+ 4r 2 + er(/(26(7 + 2n) + 7a(l + 6n)) + (17a + 46)/r) + e(76 + a(12 - 5n))(-2 + n) csc6> 2 )) 

+ 2f8 7 r 5 (er(f(-5b(2 + n) - a(14 + 25n)) - 2(9a + 56)/r) + 2(en(6 - llan - 46n) - 5r 2 

+ en(6(-5 + 4n) + a(-12 + lln)) esc 8 2 )) - A8 S r 2 (er(~2f(-8a + 6 + 19an + 46n)r 

- 4(4a + 6)/ (3) r 2 + /(a(76 + (150 - 19n)n) + 6(58 + (9 - An)n) + 6(3a + Ab)8nr 2 cot 8)) 
+ 2(en(a(44 - 15n)n + 6(-6 + 17n)) + (9 + 6n)r 2 + en(6(20 - 17n) 

+ a(50 + n(-62 + 15n))) esc 8 2 ))) + 4(-2(10a + 7b)f 3 8 Z er 6 + / 2 er 3 (2(18a + 17b)8r 

+ 8(20a(2 + n) + 66(6 + n) + (26a + 7b)8nr cot 8)) + 2n cot 8(~en(b + an) + r 2 - 2afer 2 

+ e(-2 + n)(b + an) esc 8 2 ) + 2/r(2en(2a + 6 - 2an) cot 8 
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5 ° = ~ 192(a + b^ePF ^ [ 24 ^ £ / 6 ( fl + b + a ")(' 1 + "V + ^ A 6 (6fe(a + 6 

• 3 

+ on) (-24 + n(17 + n)) + 60e(-13 + 7n)(a + 6 + on)r + 9 r 2 (4en(l + 2n)(a + 6 + art) 

+ 6(1 + n)r 2 + 3e(l + n)r(/(3(a + 6) + 5an) + 2o/r))) + 2<9 2 / 4 r 4 (-126> 3 e(a(-l + n)(-80 

+ n(-26 + 3n)) + 6(80 + n(-69 + In))) - 1440(a + b)fer 3 + 240(a + 6)06>er 2 (3fl(-3 + n) 

+ 46> (3) r) + 406>V(4en(6 + 46n + a(3 + 5ra)) + 3r 2 + 3er(3/(a + 26 + 3an) 

+ 4a/r)) + /6>V(4a£n(l + 2rc) + 6r 2 + 3er((3a + 6)/(l + n) + 2 (a + 6)/V)) 

- 246> 2 er(0(6(71 - 42n + 4n 2 ) + o(71 - n(45 + 2n))) + 4(a + 6)r(26> (3) (-3 + n) + 6> (4) r)) 
+ 46> 5 r 2 (2£7i(6(2 + lln) + a(14 - (-24 + n)n)) + 3(5 + 4n)r 2 + 3er(/(26(8 + n) 

+ a(13 + n(27 + 5n))) + r(/(5a - 36 - 2a?i) - 2a/ (3) r)))) - 2/(384(a + 6)6» (3) fer 3 

- 18(46a + 43b)6f 2 8 2 er 5 + 27(2a + 36)/ 3 6» 5 £r 7 - 2/6> 3 r 3 (2en(6(5 + 4n) 

- 3a(4 + 7n)) + 2lr 2 + 3er(/(2a(62 + n) + 6(102 + In)) + 126/r)) 

+ 40r(2en(-6 + a (4 + 7n)) - 3r 2 + 6er(/(20a(2 + n) + 6(40 + 17n)) + (10a + 96) fr)) 

+ 46(2en 2 (-b + a (4 + 7n)) - 3(2 + 7i)r 2 + 6er(/(2a(2 + n){2 + 3n) + 6(8 + n(15 + n))) 

+ r(2/(3o(2 + n) + 6(5 + n)) + (2a + 6)/ (3) r)))) + / 3 r 2 (480(a + b)(fer 3 

- 6(5a + b)6f 2 6 6 er 7 - 3(a + b)f 3 9 er 9 + 480(a + 6)^er 2 (30(4 + n) 

+ 46> (3) r) + 246» 3 e(6(48 - 13(-4 + n)n) + a(48 + n(60 + (8 - 5n)n)) - 30(a + b)0 2 fr 5 ) 

+ 2./<9 7 r 5 (4en(6 + 46?i + a(6 + lln)) + 30r 2 + 3er(./(56(2 + n) + a(14 + 25n)) 

+ +2(9a + 56) /V)) - 246> 2 er(0(-24O(a + 6) - (68a + 476)n + (13a + 166)n 2 ) 

+ 16(a + 6)r((9 (3) (-l + 2n) + 6> (4) r)) + 86>V(24(a + b)6 (3) fer 2 + 8{2en(8a + 56 + 13an 

+ 16bn) + 9r 2 + 3er(/(-14a - 236 + 43an + 16bn) + 6(3a + 26) fr))) 

- 46> 5 r 2 (4£7i(-6(l + 17n) + 5a(-5 + n(-7 + 3n))) - 18(3 + 2n)r 2 

+ 3er(/(6(-58 + n(-9 + in)) + a(-76 + ra(-150 + 19n))) + 2r(/(-8a + 6 + 19an + 46n) 

+ 2(4a + 6)/ (3) r)))) - 2/ 2 (18(7a + I0b)9 f 2 4 er 7 + 9(2a + b)f 3 6 7 er 9 

+ 246>e(n 2 (36 + a(2 + n)) - 90(a + b)6 2 fr 5 ) - 2f9 5 r 5 (2(a + Ab)en(l + 2n) + 27r 2 

+ 9er(/(56(-l + n) + a (4 + lln)) + 8(a + 6)/r)) + 24er(6(a(2 + n)(4 + 5n) 

+ 6(8 + n(15 + 4n))) + 4(a + 6)r(26» (3) (2 + n) + 6» (4) r)) + 46» 2 r 3 (72(a + b)6 {3) fer 2 

+ 0(4(a - 2b)en(l + 2n) - 9r 2 + 18er(/(a(-48 + n) + 6(-46 + 3n)) + (a + 6)/r))) 

+ 46> 3 r 2 (2£7i(-96n + a(-12 + n(-10 + 19n))) - 3(7 + 4n)r 2 + 6er(./(6(-78 + n(13 + 2«)) 

+ 2a(-42 + n(-5 + 6n))) + r(/(6(9 + 4n) + 2a(5 + In)) + (5a + 26)/ (3 V)))) 

+ 8/r(-3(18a + llb)Bfer 2 + 6(2sn{b - a (4 + 7n)) + 3r 2 + 3er(-4(2a + b)fr 

+ /(-10a(2 + 36(6 + + (10a + 76)/6> 2 r 3 )))) ] . (A.20) 
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